Crack propagation in quasicrystals 
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Crack propagation is studied in a two dimensional decagonal model quasicrystal. The simulations 
reveal the dominating role of highly coordinated atomic environments as structure intrinsic obstacles 
for both dislocation motion and crack propagation. For certain overloads, these obstacles and the 
quasiperiodic nature of the crystal result in a specific crack propagation mechanism: The crack tip 
emits a dislocation followed by a phason wall, along which the material opens up. 



Quasicrystals are characterized by an unusual struc- 
ture: They are orientationally ordered, but with noncrys- 
tallographic, e.g., five- or eightfold symmetry axes [Q. 
Their translational order is not periodic, as in conven- 
tional crystals, but quasiperiodic. In the last few years 
it has become possible to measure physical properties 
of quasicrystals, such as mass and heat transport, plas- 
ticity, and fracture, since high-quality, thermodynami- 
cally stable quasicrystals could be grown in millimeter 
sizes (|] . Because of the peculiar structure, quasicrystals 
open new possibilities for studying structure-property re- 
lations. The mechanical behavior of a solid, in particular, 
is strongly influenced by structure, because it is mainly 
governed by defects, which are related to the order. Sev- 
eral studies of the standard icosahedral quasicrystallinc 
alloys AlMnPd and AlCuFe have disclosed that the sam- 
ples are brittle and hard at room temperature S, but 
become ductile at about 80% of the melting temperature, 
and then can be deformed up to 30% without hardening 
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Concerning ductility it has been confirmed both by 
transmission electron microscopy [B and by computer 
simulations ||, that motion of dislocations is the mech- 
anism of plastic deformation in quasicrystals. The exis- 
tence of dislocations in quasicrystals is due to the fact 
that quasiperiodic structures can be considered as irra- 
tional cuts through a higher dimensional periodic crys- 
tal 0. Hence the Burgers vector of a dislocation also 
is of higher dimension, with a component in physical 
space and a component in the complementary orthogonal 
space. The dislocation is described by a "phonon"-like 
displacement field u(x) and a "phason" -like displacement 
field w(x), whose contour integrals yield 6" and b ± , re- 
spectively. 

Molecular dynamics simulations of a sheared quasicrys- 
tal have revealed that close to the melting temperature 
the phason field is spread out around the dislocation core, 
whereas at low temperatures it is confined to a plane ex- 
tending from and being bound by the dislocation line 
H|j. The phase of the quasicrystal is differing by b 1 - at 
the two sides of this wall, which therefore is denoted as 
"phason- wall" . 

In this article we demonstrate - by employing atom- 
istic simulations of fracture processes - that phason walls, 
being a consequence of the quasiperiodic structure, are a 



crucial element in the brittle fracture mechanism at low 
temperature. Furthermore, we identify atomic clusters, 
inherent to the quasicrystalline structure, as obstacles 
and scattering centers for the brittle crack. 

The study of a strongly nonlinear process like dynam- 
ical fracture in a complex material requires simplifica- 
tions. First we consider a two dimensional system (2D) 
to evade the analysis of a three dimensional crack front. 
Secondly we simulate a perfectly ordered quasicystalline 
representative of the binary tiling syste m |H |, derived 
from the Tubingen triangle tiling (TTT) (Tj. A large 
atom (L) is placed at each vertex and a small atom (S) 
into the center of each large triangle of the original tiling. 
When the centers of neighboring atoms of different type 
are connected, one obtains a rhombic tiling, which we 
denote as "bond representation" of the direct "atom rep- 
resentation" . 

Qualitatively our model shares central properties with 
structure models for real quasicrystals, which are (I) a 
decagonal diffraction pattern inherited from the TTT |ll| ] 
and (2) a hierarchical structure of clusters. These are 
built by a first shell of small atoms and a second shell 
of large atoms surrounding a central large atom (Fig. |l|), 
and can be considered as 2D versions of Mackay clusters 
that are significant for many real quasicrystals. Further- 
more, assemblies of such clusters form super clusters and 
so on. The smallest clusters are situated on five families 
of parallel lines, mutually rotated by 36°, with a small 
and a large separation within each family, arranged in a 
Fibonacci sequence (Fig. ||). 

The atoms interact with Lennard- Jones potentials 
truncated at 2.5 tq. The bond lengths are chosen to be 
the geometrical ones B . The unit of length in our sim- 
ulation ro = tls is given by the length of the LS bond. 
The bond energies are £o = £ls an d £ll = £ss = 0.5£o 
to avoid phase separation. All masses are set to 1. The 
unit of time is therefore to = ro y^m/eo. Due to their high 
coordination the decagonal clusters become the tightest 
bound structural units of the quasicrystal model for all 
choices of Sll and ess smaller than one. 

Preferred cleavage planes of the 2D model quasicrys- 
tal, i.e. the "easy lines" for fracture, are those in between 
the large separation of the parallel lines connecting the 
clusters, because a crack can follow these lines without 
cutting the tightly bound ten-rings. The surface ener- 
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gies [H of different crystallographically equivalent, but 
structurally distinct lines (1-4 in Fig. ||) are lowest for 
those between the widely spaced ten-rings (3,4 in Fig. ||), 
namely 1.9 £o/ r o as compared to 2.3 £o/^o for the other 
surfaces studied here (1,2,5 in Fig. |^). Furthermore, they 
also have the lowest "unstable stacking fault energy" , i.e. 
increase in energy when the two halfs are shifted relative 
to each other. Hence these lines are also preferred for 
dislocation motion When a dislocation is gliding 

through a quasicrystal, leaving a phason wall behind, it 
creates an even more preferred cleavage line, since the 
phason wall is characterized by a positive defect energy. 
For the shortest Burgers vector the cleavage energy along 
the phason wall can be reduced by an additional 0.3 £oAo 

& 

The fracture simulations are performed on a rectan- 
gular slab of approximately 250 000 atoms, 250 tq high 
and 1.000 rowide. The slab is uniformly strained by an 
amount A normal to the long axis, and in the outer- 
most layer of width 3ro all atoms are held fixed at these 
positions. A precrack of length 200r is inserted from 
one side along an easy line by cutting the bonds and 
relaxing the sample. For the relaxation the externally 
applied homogeneous strain is fixed at the Griffith value, 
at which the strain energy per unit length just equals 
twice the surface energy. The resulting displacement field 
then is scaled linearly above the critical dilation where 
the crack tip starts moving. The evolution of the sys- 
tem is followed by standard molecular dynamics (MD) 
technique where the atoms obey Newtonian equations of 
motion. The system was started at an initial temperature 
of T=10~ 6 T me / t . The IMD molecular dynamics package 
was used \T%. During crack propagation, the crack tip 
bonds are broken and the strain energy above the surface 
energy (the overload) is released in the form of acoustic 
emission or dislocation generation. To avoid heating of 
the slab and phonon reflections from the boundaries an 
elliptical stadium is created outside which the waves are 
damped gradually by a ramped friction term in the equa- 
tions of motion p5| . 

For the phason-free precracked sample the critical 
strain is e c =1.12%. Straight brittle cleavage along the 
easy line occurred only very close to this value and for 
short segments. In all other cases a rough surface is cre- 
ated. The propagating crack is monitored and its current 
position is plotted versus time for three different loads in 
Fig. [| Depending on the load level one can distinguish 
three different fracture regimes. 

An intermittent propagation regime is found at low 
load levels. For example, at a strain of 1.131% the crack 
proceeds with constant velocity along an easy fracture 
line until it hits an obstacle (Fig. [|, top), which can be 
a complete or an incomplete ten-ring. It then stops and 
emits a dislocation along one of the easy planes inclined 
by 36°, along which the shear stress is high. The disloca- 
tion moves away from the crack tip until it is stopped it- 
self by an obstacle(Fig. ||, middle), whereupon the entire 
configuration rests for some time. The dislocation has 



left a phason wall in its wake. This phason wall carries 
some interface energy, and hence the cohesive strength of 
the material is reduced at this wall. After some time of 
strain build-up the crack begins to open along the phason 
wall (Fig. [|, bottom). It immediately regains its initial 
propagation velocity until the next obstacle is met. Due 
to the tightly bound complete ten-rings, the cleavage sur- 
face roughness is about 10 r$. 

At intermediate loads, between 1.3% and 1.5% strain, 
a regime of steady crack propagation is observed. The 
crack tip proceeds with almost uniform velocity, but the 
appearance of the cleavage plane is unchanged. The 
dislocation-emission-phason-wall cleavage mechanism is 
still operates but without stopping the propagating 
crack. 

At even higher overloads (> 1.5% strain) we observe 
an onset of crack tip instabilities, which create defects 
that are not bound to the crack tip. At a strain of 
1.5% (Fig. [5]) we observe dislocation emission at an an- 
gle of 72°, which reduces (due to shielding) the load on 
the crack tip. Arrested by an obstacle, the dislocation 
moves back and annihilates at the fracture surface af- 
ter the crack tip has propagated further. Such a pro- 
cess, so-called "virtual dislocation emission", was first 
proposed by Brede and Haasen to explain their fracture 
experiments in silicon At a load of 1.7% the crack 
bifurcates and stops. 

In quasicrystals as in periodic crystals, we find a min- 
imal velocity for brittle crack propagation of 14% of the 
shear wave velocity vt at low loads. This suggests that 
there exists a lower band of forbidden crack velocities 
|j~7| , ^8| . With rising load, the crack tip speed increases 
to about 35% of vt- This value is again similar to the 
maximum crack tip speed in periodic crystals Jl9| . The 
lower than expected upper crack velocity, which linear 
elasticity theory would place at the Rayleigh wave veloc- 
ity |2l| [vjiayieigh 0.9 vt), may be explained with the 
non-linearity of the atomic interaction, which has been 
shown to drastically reduce crack tip velocities in peri- 
odic crystals |l^j2^]. At high loads the upper limit on 
the crack velocity is determined by the onset of crack 
tip instabilities: (virtual) dislocation emission and crack 
branching. It seems likely that the instability is driven 
by a critical overload or a critical accumulation of specific 
shear waves. 

While the crack tip velocity is of similar magnitude in 
quasicrystals and periodic crystals, the crack propaga- 
tion mode differs significantly. At small overloads an in- 
termittent crack propagation mode occurs, during which 
the average crack velocity can be low but microscopically 
separated into distinct propagation and arrest periods. 
The fracture surfaces are rough, governed by cluster in- 
homogeneities. Intermittency and roughness can both be 
attributed to the novel dislocation-emission-phason-wall 
fracture mechanism. 

At comparable external loads the brittle fracture of 
perfect periodic crystals is characterized by steady crack 
propagation along well defined flat and smooth cleav- 
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age planes [ [L8| . In comparison to periodic crystals, two 
specific structural components of the quasicrystal are im- 
portant to explain the observed behaviour. Firstly, the 
role of the ten-rings as tightly bound clusters which act 
as structure intrinsic obstacles to both dislocations and 
cracks is obvious. At low overloads the clusters may tem- 
porarily arrest the crack, and even at high overloads and 
high crack tip speeds, the clusters act as scattering cen- 
ters for the cracks. 

In perfect agreement with a previous study of dislo- 
cation mobility in quasicrystals |Dj |, we find here that 
the clusters also act as obstacles to dislocation motion. 
As a consequence, dislocations that have been nucleated 
at the crack tip cannot get farther away than the next 
sufficiently strong obstacle. Since the interaction of the 
crack tip and the dislocation produces a driving force on 
the dislocation which decreases with distance from the 
crack tip p^] , the strength of a cluster, measured as the 
capability to arrest the dislocation, increases with dis- 
tance. Strong local structural variations, which are a 
characteristic of the quasicrystalline structure, are there- 
fore very efficient in blocking dislocation motion and tend 
to keep dislocations very close to the crack tip. Usually 
this would be considered beneficial, since a dislocation 
emitted from the crack tip also shields the crack tip from 
the applied load, with the shielding decreasing with dis- 
tance to the crack tip [||| . As a consequence, one would 
expect a high fracture toughness if dislocation emission 
occurs, while dislocations are kept close to the crack tip. 
However, this is not what happens, due to the second 
structural peculiarity of quasicrystals: the phason wall 
left behind a moving dislocation. 

Since the phason wall carries excess energy, it provides 
a new low-energy path for the crack and thereby effec- 
tively allows the crack to bypass the structure-intrinsic 
obstacles. The existence of the phason wall, trailing dis- 
locations nucleated at the crack tip, is ultimately respon- 
sible for the observed brittleness of the quasicrystals. 

Although the binary model quasicrystal used in this 
study is very simplistic, it reflects the most important 
features of real quasicrystals, i.e. quasiperiodicity and 
clusters. It permits conclusions relevant to real qua- 
sicrystals. Our simulations reveal the dominating role 
of highly coordinated atomic environments as structure 
intrinsic obstacles for both crack and dislocation motion. 
Possible toughening effects, however, are over compen- 
sated by the embrittling nature of the phason wall. The 
phason wall modifies the Griffith criterion and leads to 
an immediate crack opening of the solid in the path of a 
single dislocation. The existence of the phason degree of 
freedom gives raise to this special mechanism of brittle 
fracture in quasicrystals. 
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FIG. 1. Ten-rings in atom (middle) and bond representa- 
tion (right) and their origin in the triangle tiling (left). 





FIG. 4. Crack propagation mechanism in the intermittent 
regime. Top: The crack tip has stopped. Middle: It has emit- 
ted a dislocation which is followed by a phason wall. Bottom: 
The quasicrystal has opened along the wall. 



FIG. 2. Binary tiling in the bond representation and spe- 
cial planes, whose surface energy is explained in the text. 
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FIG. 3. Propagated distance of crack tip versus time. De- 
pending on the dilation three different regimes are found: An 
intermittent (e=1.13%, 1.131% and 1.14%), smooth (s=1.3%) 
and unstable regime (e=1.5 and 1.6%). 
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FIG. 5. Virtual dislocation emission observed at a strain of 
1.5% in the bond representation. Top: Snapshot of a simula- 
tion after 73.5to; a dislocation has been emitted at an angle 
of 72°. Bottom: After 77to the dislocation has returned. 
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